CUTS AND FLOWS OF CELL COMPLEXES 
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Abstract. We study the vector spaces and integer lattices of cuts and 
flows associated with an arbitrary finite CW complex, and their rela- 
tionships to group invariants including the critical group of a complex. 
Our results extend to higher dimension the theory of cuts and flows in 
graphs, most notably the work of Bacher, de la Harpe and Nagnibeda. 
We construct explicit bases for the cut and flow spaces, interpret their 
coefficients topologically, and give sufficient conditions for them to be 
integral bases of the cut and flow lattices. Second, we determine the pre- 
cise relationships between the discriminant groups of the cut and flow 
lattices and the higher critical and cocritical groups with error terms 
corresponding to torsion (co)homology. As an application, we general- 
ize a result of Kotani and Sunada to give bounds for the complexity, 
girth, and connectivity of a complex in terms of Hermite's constant. 



1. Introduction 

This paper is about vector spaces and integer lattices of cuts and flows 
associated with a finite cell complex. Our primary motivation is the study 
of critical groups of cell complexes and related group invariants. The crit- 
ical group of a graph is a finite abelian group whose order is the number 
of spanning forests. The definition was introduced independently in several 
different settings, including arithmetic geometry [21], physics [ID] , and alge- 
braic geometry [2] (where it is also known as the Picard group or Jacobian 
group). It has received considerable recent attention for its connections to 
discrete dynamical systems, tropical geometry, and linear systems of curves; 
see, e.g., SI Q2] . 

In previous work [15] . the authors extended the definition of the critical 
group to a cell complex E of arbitrary dimension. To summarize, the critical 
group K(Ti) can be calculated using a reduced combinatorial Laplacian, 
and its order is a weighted enumeration of the cellular spanning trees of E. 
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Moreover, the action of the critical group on cellular (d — l)-cochains gives a 
model of discrete flow on £, generalizing the chip- firing and sandpile models; 
see, e.g., [HfTU]. 

Bacher, de la Harpe, and Nagnibeda first defined the lattices C and T of 
integral cuts and flows for a graph [2 J . By regarding a graph as an analogue of 
a Riemann surface, they interpreted the discriminant groups C"/C and IF* jT 
respectively as the Picard group of divisors and as the Jacobian group of 
holomorphic forms. In particular, they showed that the critical group K{G) 
is isomorphic to both C"/C and T* /F. Similar definitions and results appear 
in the work of Biggs [I] . 

In the present paper, we define the cut and flow spaces and cut and flow 
lattices of a cell complex £ by 

Cut(E) = im R d\ Flow(S) = ker R d, 

C(£) = im z «9*, J-(£)=ker z 3, 



where d and d* are the top cellular boundary and coboundary maps of £. 
In topological terms, cut- and flow-vectors are cellular coboundaries and 
cycles, respectively. Equivalently, the vectors in Cut(£) support sets of 
facets whose deletion increases the codimension-1 Betti number, and the 
vectors in Flow(S) support nontrivial rational homology classes. In the 
language of matroid theory, cuts and flows correspond to cocircuits and 
circuits, respectively, of the cellular matroid represented by the columns 
of d. 

In the higher-dimensional setting, the groups C a /C and /F are not nec- 
essarily isomorphic to each other. Their precise relationship involves several 
other groups: the critical group K(Ti), a dually defined cocritical group 
K*{Y>), and the cutflow group Z n /(C © F). We show that the critical and 
cocritical groups are respectively isomorphic to the discriminant groups of 
the cut lattice and flow lattice, and that the cutflow group mediates between 
them with an "error term" given by homology. Specifically, if dim £ = d, 
then we have the short exact sequences 

-»• Z n /(C © F) -> if (£) ^ C s /C -> T(^ d _!(S; Z)) 0, 
-»• T(# d _i(E;Z)) -4Z n /(CeJ) ^ K*(Z)^P/F ^0 



(Theorems 7.6 and 7.7) where T denotes the torsion summand. The sizes 



of these groups are then given by 



|C»/C| = |K(E)|=t(E) 
\P/F\ = \K*(E)\ =r*(E) 
\Z n /(C(BF)\ =r(E)/t 



r*(E) -t" 

r(E)/t 2 , 

r*(E)-t, 
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(Theorems 8.1 and 8.2), where t = |T(JTd_i(E;Z))| and r(S) and r*(E) are 



the weighted enumerators 

r(E) = £|T(fl d _i(T;Z))| 2 , r*(S) = £ \T(H d (Q, T; Z))| 5 



where both sums run over all cellular spanning forests T C S (see equa- 
tion p])) and is an acyclization of T (see Definition |7.3[ ). 

Before proving these results, we study the cut space (Section Hi), the flow 
space (Section [5]), and the cut and flow lattices (Section [6]) in some detail. 
In order to do this, we begin in Section [3] by describing and enumerating 
cellular spanning forests of an arbitrary cell complex, generalizing our earlier 
work \13\ Q3] . The methods and results are very close to those of Lyons [25] , 
but our technical emphasis is slightly different. For every cellular spanning 
forest TCE, the fundamental bonds of T give rise to a basis of the cut space 



(Theorem 4.7). Likewise, the fundamental circuits of T induce a basis of the 



flow space (Theorem 5.5). Moreover, under certain conditions on T, these 



bases are in fact integral bases for the cut and flow lattices (Theorems 6.1 



and 6.2) 



In the graphic case, these basis vectors are signed characteristic vectors 
of fundamental bonds and fundamental circuits; in particular, their entries 
are all and ±1. In higher dimension, the supports of basis elements are 
still fundamental bonds and fundamental circuits, but there is the additional 
subtlety that their entries may not scale to ±1, making it less clear how to 
associate a canonical cut- or flow- vector with a bond or circuit. We prove 
that the characteristic vectors which form the constructed bases of both the 
cut and flow spaces can be scaled so that their entries are torsion coefficients 



of homology groups of certain sub complexes (Theorems 4.10 and 5.3). 

The idea of studying cuts and flows of matroids goes back to Tutte [29J . 
More recently, Su and Wagner |28j define cuts and flows of a regular ma- 
troid (i.e., one represented by a totally unimodular matrix M); when M is 
the boundary matrix of a cell complex, this is the case where the torsion 
coefficients are all trivial. Su and Wagner's definitions coincide with ours; 
their focus, however, is on recovering the structure of a matroid from the 
metric data of its flow lattice. 

In the final section of the paper, we generalize a theorem of Kotani and 
Sunada [22], who observed that a classical inequality for integer lattices, 
involving Hermite's constant (see, e.g., [23]), could be applied to the flow 
lattice of a graph to give a bound for girth and complexity. We prove 



the corresponding result for cell complexes (Theorem 9.2), where "girth" 
means the size of a smallest circuit in the cellular matroid (or, topologically, 
the minimum number of facets supporting a nonzero homology class) and 
"complexity" is the torsion-weighted count of cellular spanning trees. 
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2. Preliminaries 

In this section we review the tools needed throughout the paper: cell 
complexes, cellular spanning trees and forests, integer lattices, and matroids. 

2.1. Cell complexes. We assume that the reader is familiar with the basic 
topology of cell complexes; for a general reference, see [20] , 

Throughout the paper, £ will denote a finite CW complex (which we refer 
to simply as a cell complex) of dimension d. We adopt the convention that £ 
has a unique cell of dimension —1 (as though it were an abstract simplicial 
complex); this will allow our results to specialize correctly to the case d = 1 
(i.e., that £ is a graph). We write £j for the set of i-dimensional cells in £, 
and for the i-dimensional skeleton of £, i.e., Sm = £j U £j-i U • • • U £o- 
A cell of dimension d is called a facet. 

Unless otherwise stated, every d-dimensional subcomplex r C £ will be 
assumed to have a full codimension-1 skeleton, i.e., T(d-i) = £(d_i)- Accord- 
ingly, for simplicity of notation, we will often make no distinction between 
the subcomplex Y itself and its set of facets. 

The symbol Cj(£) = Q(£; R) denotes the group of i-dimensional cellular 
chains with coefficients in a ring R. The i-dimensional cellular boundary 
and coboundary maps are respectively <9j(£; R) : Cj(£; R) — > Cj_i(£; R) and 
<9*(£; R) : Cj_i(£; R) — > C;(£; R); we will write simply di and d* whenever 
possible. 

The i" 1 reduced cellular homology and cohomology groups of S are re- 
spectively Hi(T,;R) = ker<9j/im<9i + i and H l (L;R) = ker <9* +1 /im<9*. We 
say that X is R-acyclic in codimension one if ^_i(S;i2) = 0. For a graph 
(d = 1), both Q- and Z-acyclicity in codimension one are equivalent to con- 
nectedness. The i th reduced Betti number is /3j(£) = dimi?j(S;Q), and 
the i th torsion coefficient tj(£) is the cardinality of the torsion subgroup 
T(fli(S;Z)). 

A pair of complexes r C £ induces a relative complex (£,T), with rel- 
ative homology and cohomology Hi(T,,T; R) and H l CE,T; R) and torsion 
coefficients t;(£,r) = |T(#j(£, T; Z))|. 

We will frequently use the fact that 

T(F d _ 1 (£;Z))^T(if d (£;Z)) (1) 

which is a special case of the universal coefficient theorem for cohomology 
[201 p. 205, Corollary 3.3]. 

2.2. Spanning Forests and Laplacians. Our work on cuts and flows will 
use the theory of spanning forests in arbitrary dimensions. Define a cellular 
spanning forest (CSF) of £ to be a subcomplex TCE such that T( rf _ 1 - ) = 
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S(d-i) and 

F d (T;Z) = 0, (2a) 

rank (T; Z) = rank If d _i(E;Z), and (2b) 

|T d | = |E d |-&(E) + fl,_i(E). (2c) 

These conditions generalize the definition of a spanning foreslj^of a graph G: 
respectively, it is acyclic, connected, and has n — c edges, where n and c are 
the numbers of vertices and components of G. Just as in the graphic case, 



any two of the conditions (2a), (2b), (2c) together imply the third. An 



equivalent definition is that a subcomplex T C X is a cellular spanning 
forest if and only if its ci-cells correspond to a column basis for the cellular 
boundary matrix d = f? d (E). 

In the case that E is Q-acyclic in codimension one, this definition special- 
izes to our earlier definition of a cellular spanning tree [14} Definition 2.2]. 

There are two main reasons that enumeration of spanning forests of cell 
complexes is more complicated than for graphs. First, many properties of 
graphs can be studied component by component, so that one can usually 
make the simplifying assumption of connectedness; on the other hand, a 
higher-dimensional cell complex cannot in general be decomposed into dis- 
joint pieces that are all acyclic in codimension one. Second, for complexes of 
dimension greater than or equal to two, the possibility of torsion homology 
affects enumeration. 

Define the i th up-down, down-up and total Laplacian operators^] on E by 

Lf = d i+1 d* +1 : Cj(E; R) -> d(E; R), 

L^ = d*di: C7 i (S;i?)^C i (S; J R), 
L tot = L ud + L du_ 

Moreover, define the complexity of E as 

r(E)=r i (E)= £ iTC^CrjZ))! 2 . (3) 

CSFs TCS 

The cellular matrix-tree theorem |14| Theorem 2.8] states that if S is Q- 
acyclic in codimension one and L-f is the submatrix of LH£-,(E) obtained 
by deleting the rows and columns corresponding to the facets of a (d— 1)- 
spanning tree T, then 

|T(^ 2 (E; Z ))P 
|T(ff„_ 2 (T;Z))P 



^That is, a maximal acyclic subgraph of G, not merely an acyclic subgraph containing 
all vertices. 

These are discrete versions of the Laplacian operators on differential forms of a Rie- 
mannian manifold. The interested reader is referred to [16] and [12] for their origins in 
differential geometry and, e.g., [Ill 1171 12S] for more recent appearances in combinatorics. 
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In Section|3j we will generalize this formula to arbitrary cell complexes (i.e., 
not requiring that £ be Q-acyclic in codimension one). This has previously 
been done by Lyons [22] in terms of slightly different invariants. If G is a 
connected graph, then r(G) is just the number of spanning trees, and we 
recover the classical matrix-tree theorem of Kirchhoff. 

2.3. Lattices. Starting in Section [6j we will turn our attention to lattices 
of integer cuts and flows. We review some of the general theory of integer 
lattices; see, e.g., fTj Chapter 12], pi Chapter 14], [2U Chapter IV]. 

A lattice £ is a discrete subgroup of a finite-dimensional vector space V; 
that is, it is the integer linear combinations of some basis of V. Every 
lattice £ C R n is isomorphic to Z r for some integer r < n, called the 
rank of £. The elements of £ span a vector space denoted by £ ® M. For 
£ C Z n , the saturation of £ is defined as £ = (£ ® 1) (1 Z". An integral 
basis of £ is a set of linearly independent vectors v\ , . . . , v r £ £ such that 
£ = {ci«i+- • -+c r v r : Ci G Z}. We will need the following fact about integral 
bases of lattices; the equivalences are easy consequences of the theory of free 
modules (see, e.g., [H Chapter 12], [211 Chapter IV]): 

Proposition 2.1. For any lattice £ C Z n , the following are equivalent: 

(a) Every integral basis of £ can be extended to an integral basis of Z n . 

(b) Some integral basis of £ can be extended to an integral basis of7L n . 

(c) £ is a summand of7L n , i.e., 7L n can be written as an internal direct 
sum £ CJ . 

(d) £ is the kernel of some group homomorphism Z n — > Z m . 

(e) £ is saturated, i.e., £ = £. 

(f) Z n /£ is a free 7L-module, i.e., its torsion submodule is zero. 

Fixing the standard inner product (•, •) on l n , we define the dual lattice 
of £ by 

£ , = {d££®R: {v,w)eZ Vw e £}. 

Note that £ a can be identified with the dual Z-module £* = Hom(£, Z), and 
that = £. A lattice is called integral if it is contained in its dual; for 

instance, any subgroup of Z n is an integral lattice. The discriminant group 
(or determinantal group) of an integral lattice £ is £^/£; its cardinality 
can be calculated as det M T M, for any matrix M whose columns form an 
integral basis of £. We will need the following facts about bases and duals 
of lattices. 

Proposition 2.2. |18l Section 14.6] Let M be an n x r integer matrix. 

(a) If the columns of M form an integral basis for the lattice £, then the 
columns of M(M T M)^ 1 form the corresponding dual basis for £K 

(b) The matrix P = M(M T M)~ l M T represents orthogonal projection 
from M n onto the column space of M. 

(c) If the greatest common divisor of the r x r minors of M is 1, then 
£^ is generated by the columns of P. 
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2.4. The cellular matroid. Many results of the paper may be expressed 
efficiently using the language of matroids. For a general reference on ma- 
troids, see, e.g., [27]. We will primarily consider cellular matroids. The cel- 
lular matroid of £ is the matroid .A/f (£) represented over M by the columns 
of the boundary matrix d. Thus the ground set of ,M(£) naturally cor- 
responds to the al dimensional cells X^, and .M(£) records which sets of 
columns of d are linearly independent. If £ is a graph, then .M(£) is its 
usual graphic matroid, while if £ is a simplicial complex then .M(£) is its 
simplicial matroid (see [9]). 

The bases of .A/f (£) are the collections of facets of cellular spanning forests 
of £. If r is the rank function of the matroid .M(£), then for each set of facets 
B C £ d , we have r{B) = ranked, where 8b is the submatrix consisting of 
the columns indexed by the facets in B. Moreover, we have 

r(£) := r(£ d ) = rankTW(S) = rank d = \E d \ - A*(£) 

by the definition of Betti number. 

A set of facets B C £^ is called a cut if deleting B from £ increases its 
codimension-one homology, i.e., /?d__i(£ \ B) > /3^„i(£). A cut B is a frond 
if r(£ \ B) = r(£) - 1, but r((£ \ S) Ua) = r(£) for every a £ B. That is, 
a bond is a minimal cut. In matroid terminology, a bond of £ is precisely 
a cocircuit of .M(£), i.e., a minimal set that meets every basis of .A4(£). 
Equivalently, a bond is the complement of a flat of rank r(£) — 1. If T is a 
cellular spanning forest (i.e., a basis of .M(£)) and a G is a facet, then 
the fundamental bond of the pair (T, a) is 

bo(Tf» = o-U{pG£ d \T:T\o-U / oisa CSF} . (4) 

This is the fundamental cocircuit of the pair (T,<r) of .M(£) [271 p. 78]. 

While we will frequently adopt a matroid-theoretic perspective, it is im- 
portant to point out that the cut and flow spaces and lattices of a complex £ 
are not matroid invariants, i.e., they are not determined by the cellular ma- 
troid Ai(£). (See [28] for more on this subject.) Below is a table collecting 
some of the standard terminology from linear algebra, graph theory, and 
matroid theory, along with the analogous concepts that we will be using for 
cell complexes. 



Linear algebra 


Graph 


Matroid 


Cell complex 


Column vectors 


Edges 


Ground set 


Facets 


Independent set 


Acyclic subgraph 


Independent set 


Acyclic subcomplex 


Min linear dependence 


Cycle 


Circuit 


Circuit 


Basis 


Spanning forest 


Basis 


CSF 


Set meeting all bases 


Disconnecting set 


Codependent set 


Cut 


Min set meeting all bases 


Bond 


Cocircuit 


Bond 


Rank 


# edges in spanning forest 


Rank 


# facets in CSF 



Here "codependent" means dependent in the dual matroid. 
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3. Enumerating Cellular Spanning Forests 

In this section, we study the enumerative properties of cellular spanning 
forests of an arbitrary cell complex S. Our setup is essentially the same as 
that of Lyons |25t §6], but the combinatorial formulas we will need later, 



namely Propositions 3.2 and|3.4[ are somewhat different. As a corollary, we 



obtain an enumerative result, Proposition 3.5, which generalizes the simpli- 
cial and cellular matrix-tree theorems of [13] and p3| (in which we required 
that £ be Q-acyclic in codimension one). The result is closely related, but 
not quite equivalent, to Lyons' generalization of the cellular matrix-tree the- 
orem \25\ Corollary 6.2]. 

The arguments require some tools from homological algebra, in particular 
the long exact sequence for relative homology and some facts about the 
torsion-subgroup functor. The details of the proofs are not necessary to 
understand the constructions of cut and flow spaces in the later sections. 

Let £ be a d- dimensional cell complex with rank r. Let T C £ be a sub- 
complex of dimension less than or equal to d — 1 such that T(d-2) = ^(d-2)- 
Thus the inclusion map i: T — > £ induces isomorphisms i* : Hk(T;Q) — > 
H k (Z;Q) for all k < d- 2. 

Definition 3.1. The subcomplex r C £ is called relatively acyclic if in fact 
the inclusion map i : T — > £ induces isomorphisms : Hk(T; Q) — > -#&(£; Q) 
for all k < d. 



By the long exact sequence for relative homology, T is relatively acyclic if 
and only if H d (T,; Q) — > H d (E, T; Q) is an isomorphism and -fffc(£, T; Q) = 
for all k < d. These conditions can occur only if |r^_i| = — r. This 

quantity may be zero (in which case the only relatively acyclic subcomplex 
is S( rf _ 2 )). A relatively acyclic subcomplex is precisely the complement of a 
(d — l)-cobase (a basis of the matroid represented over R by the rows of the 
boundary matrix d) in the terminology of Lyons [25j. 

Two special cases are worth noting. First, if d = 1, then a relatively 
acyclic complex consists of one vertex in each connected component. Second, 
if Hd-i(T,;Q) = 0, then T is relatively acyclic if and only if it is a cellular 
spanning forest of 

For a matrix M, we write Ma,b for the restriction of M to rows indexed 
by A and columns indexed by B. 

Proposition 3.2. Let T C T C S be subcomplexes such that dimT = d; 
dimT = d-l; \T d \ = r; \T d -i\ = ~r; ^(d-i) = andT( d _ 2 ) = 

£(d-2)- Also, let R = Yid-i \ r. Then the following are equivalent: 

(a) The r x r square matrix d = c^y is nonsingular. 

(b) # d (T,r ; Q) = and # d _i(T, r ; Q) = 0. 

(c) H d (T,T;Q) = or H d -i(T,T;Q) = 0. 

(d) T is a cellular spanning forest of T, and T is relatively acyclic. 
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Proof. The cellular chain complex of the relative complex (T, T) is 



o^c d (T,rv 



o 



with other terms zero. If d is nonsingular, then H d (T, T; Q) and H d _\ (T, T; <Q] 
are both zero; otherwise, both are nonzero. This proves the equivalence of 
(a), (b) and (c). 

Next, note that H d (T; Q) = (because T has no cells in dimension d) and 
that H d _ 2 (T, T; Q) = (because T( d -2) = ^(d-2))- Accordingly, the long 
exact sequence for relative homology of (T,T) is 







H d (T;Q 
H d -i(F; 



-> #d-i(T; 
-> ^- 2 (T; 



#d-i(T,r ; < 
o. 



(5) 



If ff d (T, r; Q) = H d -i(T, T; Q) = 0, then # d (T; Q) = (which says that T 
is a cellular spanning forest) and the rest of §5§ splits into two isomorphisms 
that assert precisely that T is relatively acyclic (recall that H d _i(T;Q) = 
H d -i(Ti;Q) when T is a cellular spanning forest). This implication is re- 
versible, completing the proof. □ 

The torsion subgroup of a finitely generated abelian group A is defined as 
the subgroup 

T(A) = {x G A: kx = for some k £ Z}. 

Note that A = T(A) if and only if A is finite. The torsion functor T is left- 
exact [21, p. 179]. Moreover, if^^^^C^Ois exact and A = T(A), 
then T(A) -)• T(B) -)• T(C) ^ is exact. We will need the following 
additional fact about the torsion functor. 

Lemma 3.3. Suppose we have a commutative diagram of finitely generated 
abelian groups 



(6) 



such that both rows are exact; A, A' are free; a is an isomorphism; (3 is sur- 
jective; and C, C are finite. Then there is an induced commutative diagram 








TB © G — TC 



• TD - 

I 

TD' 








(7) 



> TB' © G TC' > TD' > TE' 
such that G is finite and both rows are exact. Consequently 

\TB\ ■ \TC'\ ■ \TD\ ■ \TE'\ = \TB'\ ■ \TC\ ■ \TD'\ ■ \TE\ 



(8) 
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Proof. First, let Z be a maximal-rank free summand of D, i.e., D = Z © 
TD. Then im(/i) C TD. By exactness at E, the map j must map Z 
isomorphically onto a free summand of E. Hence ker j C TD as well, and 
replacing D, E with their torsion summands preserves exactness. The same 
argument implies that we can replace D' , E' with TD' , Ti?'. 

Second, note that A, A' , B, B' all have the same rank (since the rows are 
exact, C, C are finite, and a is an isomorphism). Hence f(A) is a maximal- 
rank free submodule of B; we can write B = TB © i 7 , where F is a free 
summand of B containing f(A). Likewise, write B' = TB' ®F' , where F' is 
a free summand of B' containing f'(A'). Meanwhile, /3 is surjective, hence 
must restrict to an isomorphism F —> F', which induces an isomorphism 
F/f(A) —7- F'/f'(A'). Abbreviating this last group by G, we obtain the 
desired diagram Q . Since ker g = im / C F, the map g : TB © G — > TC is 
injective, proving exactness of the first row; the second row is exact by the 
same argument. Exactness of each row implies that the alternating product 
of the cardinalities of the groups is 1, from which the formula Q follows. □ 

Proposition 3.4. Let T, be a d- dimensional cell complex, let T C £ be 
a cellular spanning forest, and let T C S be a relatively acyclic (d — 1)- 
subcomplex. Then 

td-i(T) td_i(E, r) = td_!(s) w(r, r). 

Proof. The inclusion TCE induces a commutative diagram 
^ H d -lF; Z) H d _i(Y; Z) -£ iT^T, T; Z) > # d - 2 (r ; Z) F d _ 2 (T; Z) 



V 

o ^ ff d -!(r ; z) -i fl-d_i(S;Z) 4: ^(E.r; z) > # d _ 2 (r ; z) > # d _ 2 (s ; z) 

whose rows come from the long exact sequences for relative homology. (For 
the top row, the group Hd(T ,T;Z) is free because dimT = d, and on the 
other hand is purely torsion by Proposition |3.2[ so it must be zero. For 
the bottom row, the condition that T is relatively acyclic implies that 
is an isomorphism over Q; therefore, it is one-to-one over Z.) The groups 
i/(i_i(T, T; Z) and _£?d_i(£, T; Z) are purely torsion. The first, fourth and 
fifth vertical maps are isomorphisms (the last because T( rf _ X ) = and 
the second is a surjection by the relative homology sequence of the pair 
(X, Y) (since the relative complex has no cells in dimension d — 1.) The 



result now follows by applying Lemma 3.3 and canceling like terms. □ 



As a consequence, we obtain a version of the cellular matrix-forest theo- 
rem that applies to all cell complexes (not only those that are Q-acyclic in 
codimension one). 

Proposition 3.5. Let £ be a d-dimensional cell complex and let T C S 
be a relatively acyclic {d — 1) -dimensional subcomplex, and let Lr be the 
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r d (E) = - d -;; ;, 2 detL r . 



restriction o/LJJ_ 1 (S) to the (d — l)-cells ofT. Then 

t rf -i(S) 2 

tet-i(s,r) 

Proof. By the Binet-Cauchy formula and Propositions |3.2| and |3.4[ we have 
det L r = det dr&r = Yl (deta r , T ) 2 = ^ *d-i(T,r) 2 

TCE d : |T|=r(S) CSFs TCE d 

t d _i(s,r) 2 ^ 2 t^i(s,r) 2 

= t d _!(S) 2 ^ td " l(T) = t d _!(E) 2 Td(E) 

v ; CSFs TCS d o J-V y 

and solving for Td(E) gives the desired formula. □ 

If i?d_i(E;Z) = T(i^-i(E; Z)), then the relative homology sequence of 
the pair (E, T) gives rise to the exact sequence 

-> T(^_!(E; Z)) -> T(^_!(E, T; Z)) -> T(H d _ 2 (T; Z)) -> T(if d _ 2 (S; Z)) 



which implies that trf_ 1 (E)/trf_ 1 (E, r) = t ( i_ 2 (E)/trf_ 2 (r)] so Proposition 3.5 



becomes the formula Trf(E) = t ~ 2 2 ^2 detLr- This was one of the original 
versions of the cellular matrix-tree theorem |14[ Theorem 2.8(2)]. 

Remark 3.6. Lyons |25| Corollary 6.2] proves a similar matrix- forest the- 
orem in terms of an invariant t' defined below. He shows that each row 
of ^ induces the corresponding row of ([7]). This does not quite imply 



Lemma 3.3, since one still needs to identify the "error terms" G in the 



top and bottom rows of ([7]). Doing so would amount to showing that 



3.4 



Hd-i(Y)/ ker^*) = Hd-i(E)/ker(j*) in Alternatively, Proposition 
would follow from |25|, Lemma 6.1] together with the equation 

td-aCrjt^jCrytd-aCE) = t d _ 1 (s,r)/t d _ 1 (E) 

where f = X {d _ x) \ Y and t^f) = | ker<9 d _i(E;Z)/((ker<9 d _i(E;Z) n 
ima d (£;Q)) + kercVi(r;Z))|. 

4. The Cut Space 

Throughout this section, let E be a cell complex of dimension d and rank r 
(that is, every cellular spanning forest of E has r facets). For each i < d, 
the i-cut space and i-flow space of E are defined respectively as the spaces 
of cellular coboundaries and cellular cycles: 

Cuti(E) = im(<9*: CVi(E,M) Cj(E,R)), 

Flowi(E) = ker(<9 i: Ci(E,R) -> Cj_i(E,R)). 

We will primarily be concerned with the case i = d. For i = 1, these are the 
standard graph-theoretic cut and flow spaces of the 1-skeleton of E. 

There are two natural ways to construct bases of the cut space of a graph, 
in which the basis elements correspond to either (a) vertex stars or (b) the 
fundamental circuits of a spanning forest (see, e.g. [181 Chapter 14]). The 
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former is easy to generalize to cell complexes, but the latter takes more 
work. 

First, if G is a graph on vertex set V and R is a set of ("root") vertices, 
one in each connected component, then the rows of d corresponding to the 
vertices V \ R form a basis for Cuti(G). This observation generalizes easily 
to cell complexes: 

Proposition 4.1. A set of r rows of d forms a row basis if and only if the 
corresponding set of (d — 1)- cells is the complement of a relatively acyclic 
{d — l)-subcomplex. 



This is immediate from Proposition 3.2 Recall that if Hd-i(T,;Q) = 0, 



then "relatively acyclic (d — l)-subcomplex" is synonymous with "spanning 



tree of the {d — l)-skeleton" . In this case, Proposition 4.1 is also a conse- 
quence of the fact that the matroid represented by the rows of dd is dual to 
the matroid represented by the columns of dd-i [14, Proposition 6.1]. 

The second way to construct a basis of the cut space of a graph is to fix a 
spanning tree and take the signed characteristic vectors of its fundamental 
bonds. In the cellular setting, it is not hard to show that each bond supports 
a unique (up to scaling) vector in the cut space (Lemma |4.3| ) and that 
the fundamental bonds of a fixed cellular spanning forest give rise to a 
vector space basis (Theorem 4.7). (Recall from Section [2] that a bond in 



a cell complex is a minimal collection of facets whose removal increases 
the codimension-one homology, or equivalently a cocircuit of the cellular 
matroid.) The hard part is to identify the entries of these cut- vectors. For 
a graph, these entries are all or dbl. In higher dimension, this need not 
be the case, but the entries can be interpreted as the torsion coefficients of 
certain subcomplexes (Theorem 4.10). In Section [5j we will prove analogous 



results for the flow space. 

4.1. A basis of cut-vectors. Recall that the support of a vector v = 
(yi,..., v n ) £ R n is the set 

supp(u) = {i G [n] : Vi ^ 0}. 

Proposition 4.2. \27\ Proposition 9.2.4]. Let M be a r x n matrix with 
rowspace V C W 1 , and let A4 be the matroid represented by the columns 
of M . Then the cocircuits of Ai are the inclusion-minimal elements of the 
family Supp(y) := {supp(u) : v £ V \ {0}}. 

Lemma 4.3. Let B be a bond o/S. Then the set 

Cut B (E) = {0} U{d£ Cut d (E) : supp(w) = B} 

is a one- dimensional subspace o/Cut^(S). That is, up to scalar multiple, 
there is a unique cut-vector whose support is exactly B. 

Proof. Suppose that v, w are vectors in the cut space, both supported on B, 
that are not scalar multiples of each other. Then there is a linear combina- 



tion of v, w with strictly smaller support; this contradicts Proposition 4.2 
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On the other hand, Proposition 4.2 also implies that Cuts(E) is not the 



zero space; therefore, it has dimension 1. □ 

We now know that for every bond B, there is a cut- vector supported on B 
that is uniquely determined up to a scalar multiple. As we will see, there 
is a choice of scale so that the coefficients of this cut-vector are given by 



certain minors of the down-up Laplacian L = L^ U (S) = d*d (Lemma 4.5); 



these minors (up to sign) can be interpreted as the cardinalities of torsion 



homology groups (Theorem 4.10). 



In choosing a scale, the first step is to realize the elements of Cuts(E) 
explicitly as images of the map d* . Fix an inner product (•, •) on each chain 
group Ci(Y,;M) by declaring the i-dimensional cells to be an orthonormal 
basis. (This amounts to identifying each cell with the cochain that is its 
characteristic function.) Thus, for a E Cj(£;]R), we have supp(a) = {a E 
Sj : (a, a) / 0}. Moreover, for all (3 E Cj_i(£;R), we have by basic linear 
algebra 

(da,(3) = (a,d*(3). (9) 

Lemma 4.4. Let B be a bond of E and let U be the space spanned by 
{da: <r £ B}. In particular, U is a subspace o/imfl of codimension 
one. Let V be the orthogonal complement ofU in imd, and let v be a nonzero 
element ofV. Then supp((9*v) = B. 

Proof. First, we show that 9*« / 0. To see this, observe that the column 
space of d is U + Mv, so the column space of d*d is d*U + Md*v. However, 
rank((9*<9) = rank 3 = r, and dimU = r — 1; therefore, d*v cannot be the 
zero vector. Second, if a G T,d\B, then da E U, so (d*v , a) = (v, da) = 0. It 



follows that supp(9*t>) C B, and in fact supp(<9*u) = B by Proposition 4.2 



□ 

Given a bond B, let A = {ai, . . . , a r -i} be a cellular spanning forest 
of £ \ B. Fix a facet a = a r G B, so that A U a is a cellular spanning forest 
of E. Define a vector 

r 

v = v A , a = Y,(-iy(detLf tAUuV] )da 3 G C^frZ) 

3=1 



so that 



d*v = ^(-l) J (detLj AUaVj )L d V, G Cut d (S). (10) 

3=1 



Lemma 4.5. For the cut-vector d*v defined in equation (10), 



d*v = ^2(detL% pAUa )p. 

peB 



In particular, supp((9*-u) = B. 
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Proof. For each p £ B, 



3=1 

= j2(-l) i detL% AUaV .(da j ,dp) 

3=1 

= ^(-lV'detL^ U(TVj L^ 
3=1 

= det -C/Aup.Aucr' 

where the last equality comes from expanding the row corresponding to p. 
Note that detL% p AUa / for p = a, so d*v ^ 0. On the other hand, 
by Cramer's rule, v is orthogonal to da\, . . . , <9cr r -i, so in fact (d*v, p) = 
for all p €z Y^d \ B. This establishes the desired formula for d*v, and then 



supp(9*f ) = B by Lemma 4.4 □ 



Equation (10) does not provide a canonical cut-vector associated to a 
given bond B, because d*v depends on the choice of A and a. On the other 
hand, the bond B can always be expressed as a fundamental bond bo(T, a) 
(equivalently, fundamental cocircuit; see equation Q in Section 2.4) by 
taking a to be an arbitrary facet of B and taking T = A U a, where A 
is a maximal acyclic subset of E \ B. This observation suggests that the 
underlying combinatorial data that gives rise to a cut-vector is really the 
pair (T, a). 

Definition 4.6. Let T = {o~i, 02, . . . , o~ r } be a cellular spanning forest of E, 
and let a = &i G T. The (uncalibrated) characteristic vector of the bond 
bo(T, a) is: 

x(T,a) = ^(-^(dfitL^^)^ 
3=1 



By Lemma 4.5 taking A = T \ a, we have 

pebo(T,a) 

a cut-vector supported on bo(T, a). 

The next result is the cellular analogue of |18| Lemma 14.1.3]. 

Theorem 4.7. The family {x(T,<r): a G T} is a ~SL-vector space basis for 
the cut space of E . 

Proof. Let a E T. Then suppxC^j cr) = bo(Y,c) contains cr, but no other 
facet of T. Therefore, the set of characteristic vectors is linearly indepen- 
dent, and its cardinality is | = r = dimCutd(E). □ 



CUTS AND FLOWS OF CELL COMPLEXES 



15 





(b) 




Figure 1. (a) The bipyramid 0. (b) A simplicial spanning 
tree (unfolded), (c) Deleting the bond {123, 134, 125}. 



Example 4.8. The equatorial bipyramid is the two-dimensional simplicial 
complex 9 with facet set {123,124,125,134,135,234,235} (Figure Qa)). 
Let T be the simplicial spanning tree with facets {123,124,234,135,235} 
(unfolded in Figure [T|b)). Then 

bo(T, 123) = {123, 125, 134}, bo(T, 124) = {124, 134}, 

bo(T, 135) = {135, 125}, bo(T, 234) = {234, 134}, 

bo(T,235) = {235,125}. 

In each case, the removal of the bond leaves a 1-dimensional hole (as shown 
for the bond {123, 125, 134} in Figure [ijc) ) . By Theorem [47, we have 

X(T, 123) = 75([123] + [125] - [134]), X (T, 124) = 75([124] + [134]), 
X(T, 135) = 75Q125] + [135]), X (T, 234) = 75([134] + [234]), 

X(T,235) = 75([235]-[125]), 
which indeed form a basis for Cut2(©). 

4.2. Calibrating the characteristic vector of a bond. The term "char- 
acteristic vector" suggests that the coefficients of x(T, a) should all be 



or 1, but this is not necessarily possible, even by scaling, as Example 4.12 



below will show. We would like to define the characteristic vector of a bond 
so that it carries combinatorial or topological information, avoiding extra 
factors such as the 75 in Example |4.8| We will show that the number 

td_i(x,r) 2 



t d -i(X) 2 



11 



where the sum runs over all relatively acyclic (d — l)-subcomplexes r C S, 
divides every coefficient of the characteristic vector. Moreover, we will show 
that the entries of t^tx(T, a) are (up to sign) the torsion coefficients of the 
cellular forests {T \ a U p) for p £ bo(T, a). 
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Let e^ a , be the relative sign of da, da' with respect to dA; that is, it is 
+1 or —1 according to whether da and da 1 lie on the same or the opposite 
sides of the hyperplane in im d spanned by dA. In the language of oriented 
matroids, this sign is simply a product of the entries corresponding to da and 
da 1 in one of the cocircuits corresponding to the hyperplane and determines 
the relative signs of a basis orientation on A U a and A U a' [61 Section 3.5]. 



Proposition 4.9. Let J] be a cell complex of rank r. Let T = A U a and 

T' = A U a 1 be d- dimensional cellular spanning forests of E. Let L T U T , be 
the restriction of the down-up Laplacian L dn = d*d to the rows indexed by 
T and the columns indexed by T'. Then 

detL r n r , = det<9 x cV = e^ T t d _i(T') 

where e = , 



(7,(7 



Proof. By the Binet-Cauchy formula, we have 

det L r y, = ]T(det d Tj5 )(det d s ,y>) = £(det d s ,r)(det 9 S ,T')> (12) 
s s 
where the sum runs over all sets S C Erf_i with \ S\ = r, and ds,r is the corre- 
sponding r x r submatrix of 5. By Proposition 3.2 , det c^x is nonzero if and 
only if = \ S is relatively acyclic, and for those summands Propo- 
sition IO implies |detds,x| = t d _i(T,r s ) = t (2 _i(T)t d _ 1 (E, rs)/trf_i(E). 



The sign of each summand in equation ([12j) is just e, so 
detL T u T , = y(deta 5 ,x)(deta SiT ') 



s 

= trf - l(s) A w(e) 

j . t rf _i(T)t d _i(T / ) \ - 2 
= £_ t^(E)^V 

= ^ T t d -i(T / ). □ 

Theorem 4.10. Let B be a bond. Fix a facet a £ B and a cellular forest 
A C Erf \ B , so that in fact B = bo{A U a, a). Define the characteristic 
vector of B with respect to A as 

Xa(B) := — X (AUa,a) = V e* p t d -M U p) p. 
fly ,H 

Then xa(B) is in the cut space o/E, and has integer coefficients. Moreover, 
it depends on the choice of a only up to sign. 



Proof. Apply the formula of Proposition 4.9 to the formula of Definition 4.6 



for the characteristic vector, and factor the integer py out of every coeffi- 
cient. Meanwhile, replacing a with a different facet a' G B merely multiplies 
all coefficients by e^,Je^ p = e^ a , E {±1}. □ 
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If T is a cellular spanning forest of £ and a E T4, then we define the 
characteristic vector of the pair (T, a) by 

X(T£» := X T\a(bo(T,a)) = — x(T£». (13) 

Example 4.11. Let be the bipyramid of Example |4.8| Every cellular 
spanning forest T C © is torsion-free. Moreover, the relatively acyclic sub- 



complexes r that appear in equation (11) are the spanning trees of the 1- 
skeleton 0^ (see Section |3j), which is the graph K5 with one edge removed; 
Accordingly, we have /^x = r (@(i)) = 75, so the calibrated characteristic 
vectors are as in Example |4.8[ with all factors of 75 removed. 

On the other hand, fly is not necessarily the greatest common factor of the 
entries of each uncalibrated characteristic vector, as the following example 
illustrates. 

Example 4.12. Consider the cell complex £ with a single vertex v, two 
1-cells e\ and e 2 attached at v, and four 2-cells attached via the boundary 
matrix 

02 03 05 07 

ei / 2 3 



d - 

e 2 V 5 7 

Let B be the bond {0-2,03}, so that the obvious candidate for a cut-vector 
supported on B is the row vector [2 3 0] . On the other hand, taking 
A = {05} (a cellular spanning forest of £ \ B), the calibrated characteristic 
vector of Theorem 14.101 is 

Xa(B) = [10 15 0] . 



is 



For T = AL) {02}, the uncalibrated characteristic vector of Definition 4.6 

x(T,a 2 ) = [100 150 0] . 

On the other hand, the calibration factor fiy is not gcd(100, 150) = 50, but 
rather 10, since ti(T) = 10 and the summation of equation ( |11[ ) has only 
one term, namely T = £(o)- Similarly, for A' = {07} and T' = A' U {02}, we 
have 

XA'(B)=[U 21 0], x(T', 0-2) = [196 294 0] , fi T ' = 14. 

Remark 4.13. Torsion plays a role in the characteristic vectors of bonds, 
even when S is a simplicial complex. For example, let £ be the the complete 
2-dimensional simplicial complex on 6 vertices and let T be the triangulation 
of MP 2 obtained by identifying opposite faces in an icosahedron. Then T is 
a cellular spanning forest of £ (and in fact £ has twelve spanning forests of 
this kind). For any facet a 6 T, we have bo(T, a) = £2 \ ^2 U {o~}, and the 
entries of the calibrated cut-vector include both ±2 (in position a) as well 
as rtl's (in positions £ \ T). 
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Remark 4.14. When S is a graph and T is a spanning forest, is just 
the number of vertices of S. Then, for any edge a in T, the vector xt(c) is 
the usual characteristic vector of the fundamental bond bo(T, a). 



Remark 4.15. Taking T = T in the calculation of Proposition 4.9 gives 
the equality 

E^-i(^ r ) 2 = ^||JdetL d T " 



which can be viewed as a dual form of Proposition 3.5 , enumerating relatively 



acyclic (d — l)-subcomplexes, rather than cellular spanning forests. 

5. The Flow Space 

In this section we describe the flow space of a cell complex. We begin by 
observing that the cut and flow spaces are orthogonal to each other. 

Proposition 5.1. The cut and flow spaces are orthogonal complements un- 
der the standard inner product on Crf(S;R). 

Proof. First, we show that the cut and flow spaces are orthogonal. Let 
a G Cutd = im<9^ and f3 £ Flow^ = kerc^. Then a = <9*7 for some (d — 1)- 
chain 7, and (a, /3) = (d*j,(3) = (7,9/3) = by equation ([9]). 

It remains to show that Cut^ and Flow^ have complementary dimen- 
sions. Indeed, let n = dimCd(X;]R); then dim Flow ^ = dim ker <9^ = n — 
dim \mdd = n — dim im <9*j = n — dim Cut^. □ 

Next we construct a basis of the flow space whose elements correspond 
to fundamental circuits of a given cellular spanning forest. Although cuts 
and flows are in some sense dual constructions, it is easier in this case to 



work with kernels than images, essentially because of Proposition 2.1 As a 
consequence, we can much more directly obtain a characteristic flow vector 
whose coefficients carry topological meaning. 

We need one preliminary result from linear algebra. 

Proposition 5.2. Let N be a matrix with c columns and r rows such that 
rank A = c — 1 . In particular, r > c — 1 and dim ker N = 1 . Then ker N 
has a spanning vector v = (v\, . . . , v c ) such that 

Vi = ±|T(coker Ni)\ 

where Ni denotes the submatrix of N obtained by deleting the i th column. 
In particular, Vi 7^ for all i. 

Proof. Let Q be a r x r matrix whose first r — (c — 1) rows form a Z- module 
basis for ker(A T ), and whose remaining c — 1 rows extend it to a basis 



of Z r (see Proposition 2.1). Then Q is invertible over Z, and the matrix 
P = QN = {N T Q T ) T has the form 



CUTS AND FLOWS OF CELL COMPLEXES 



19 



where M is a (c — 1) x c matrix whose column matroid is the same as 
that of N. Then keriV = kerP = kerM. Meanwhile, by Cramer's rule, 
kerM is the one-dimensional space spanned by v = (v%, . . . , v c ), where 
Vi = (— l)*detMf = ±|cokerMj| = ±|T(coker Pj)|. For each i, we have 
a commutative diagram 







im Nj 







imPr 



im.Ni 



Y 

im P~ 



T(cokeriV T ) 



T(cokerP*; 











where the hat denotes saturation (see Section 2.3). Since Q is invertible, 
it induces isomorphisms as labeled, so by the snake lemma, the third ver- 
tical arrow is an isomorphism as well. Hence |T(cokerPj)| = |T(coker A^)|, 
completing the proof. □ 

Recall that a set of facets C C is a circuit of the cellular matroid 
M(Ti) if and only if it corresponds to a minimal linearly dependent set of 
columns of 84. Applying Proposition 5.2 with N = dc (i.e., the restriction 



of d to the columns indexed by C), we obtain a flow vector whose support 
is exactly C. We call this the characteristic vector <p(C). 

Theorem 5.3. Let C be a circuit of the cellular matroid At(S), and let 
A C S be the subcomplex S(d_i) U C. Then 

V (C) = £±t d _i(A\(7)<r. 

o-ec* 

Proof. For a G C, let denote iV with the column a removed. By Propo- 

C d _i(£;Z) 



sition 5.2, it suffices to show that the two groups 

ker<9 d _i 



H d _ 1 (A\a;: 



coker Na- 



irn No- im No- 

have the same torsion summands. But this is immediate because kerc^_i is 
a summand of Cd_i(X;Z) as a free Z-module. □ 

Example 5.4. Consider the cell complex S with two vertices v\ and V2, 
three one-cells e\, e%, and es, each one with endpoints v\ and V2, and three 
two-cells ax, 0-2, and o~3 attached to the 2-cells so that the 2-dimensional 
boundary matrix is 



8 



The only circuit in S is the set C of all three 2-cells; then tp(C) = 2o~\ + 2o~2 + 
4(73, because ^ 2 (A \ <n; Z) = P^ 2 (A \ a 2 ; Z) = Z © Z 2 , but ^ 2 (A \ ct 3 ; Z) = 

z©z 2 ©z 2 . 





o"i 


2 


0"3 


61 1 


^2 


2 





e 2 







1 


e3 


V o 2 


-2 


-1 
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For a cellular spanning forest T and facet a Y, let ci(T, a) denote the 
fundamental circuit of a with respect to Y, that is, the unique circuit in 
TUd. 

Theorem 5.5. Let E be a cell complex and T C E a cellular spanning 
forest. Then the set 

Mci(T,<x)):<r0T} 
forms a M-vector space basis for the flow space of E. 

Proof. The flow space is the kernel of a matrix with |E^| columns and rank 
| Yd |, so its dimension is |E^| — |T<j|. Therefore, it is enough to show that the 
(p(c\(Y, er)) are linearly independent. Indeed, consider the matrix W whose 
rows are the vectors y>(ci(T, <r)); its maximal square submatrix W' whose 
columns correspond to E \ T has nonzero entries on the diagonal but zeroes 
elsewhere. □ 



Example 5.6. Recall the bipyramid of Example |4.8[ and its spanning tree 
T. Then ci(T, 125) = {125, 123, 135, 235}, and ci(T, 134) = {134, 123, 124, 234}. 
If we instead consider the spanning tree T' = {124, 125, 134, 135, 235}, then 
ci(T', 123) = {123, 125, 135, 235}, and ci(T', 234) = {234, 124, 125, 134, 135, 235}. 
Each of these circuits is homeomorphic to a 2-sphere, and the corresponding 
flow vectors are the homology classes they determine. Furthermore, each 
of {ci(T, 125), ci(T, 134)} and {ci(T', 123), ci(T', 234)} is a basis of the flow 
space. 

6. Integral Bases for the Cut and Flow Lattices 
Recall that the cut lattice and flow lattice of E are defined as 
C = C(E) = im z d* d C Z n , F = J"(E) = ker z d d C Z n '. 
In this section, we study the conditions under which the vector space bases 



of Theorems 4.10 and 5.5 are integral bases for the cut and flow lattices 
respectively. 

Theorem 6.1. Suppose that E has a cellular spanning forest T such that 
Hd-i(T;7j) is torsion-free. Then 

{x(T,a):aer} 

is an integral basis for the cut lattice C(E), where x(Y,cr) is defined as in 



equation (13). 

Proof. Consider the n x r matrix with columns x(Y,<r) for a G Y^. Its 
restriction to the rows is diagonal, and by Theorem |4. 7| and the hypoth- 
esis on Hd— i(T;Z), its entries are all ±1. Therefore, the x(Y,cr) form an 
integral basis for the lattice Cutrf(E) n Z n . Meanwhile, 

(Cut d (E)nZ n )/C rf (E) = T(^ d (E;Z)) T{H d ^(E; Z)) 

where the first equality is because Cutd(E) n Z n is a summand of Z n , and 
the second one is equation Q. On the other hand, ^_i(E; Z) is a quotient 



CUTS AND FLOWS OF CELL COMPLEXES 



21 



of Hd-i(T;Z) of equal rank; in particular, T(/^_i(X; Z)) = and in fact 

Cut d (s)nz n = c d (s). □ 

Next we consider integral bases of the flow lattice. For a circuit C, define 

<p(C) = -cp(C) 
9 

where <p(C) is the characteristic vector defined in Section [5] and g is the gcd 
of its coefficients. Thus <p(C) generates the rank-1 free Z-module of flow 
vectors supported on C. 

Theorem 6.2. Suppose that S has a cellular spanning forest T such that 
Hd-i(T;7*) = Hd-i(Ti;X). Then {<£>(ci(T, a)) : a T} is an integral basis 
for the flow lattice ,F(£). 

Proof. By the hypothesis on T, the columns of d indexed by the facets 
in T form a Z-basis for the column space. That is, for every a T, the 
column d a is a Z- linear combination of the columns of T; equivalently, there 
is an element w a of the flow lattice, with support ci(T, er), whose coefficient 
in the a position is ±1. But then w a and <£(ci(Y, er)) are integer vectors with 
the same linear span, both of which have the gcd of their entries equal to 1; 
therefore, they must be equal up to sign. Therefore, retaining the notation 



of Theorem 5.5, the matrix W' is in fact the identity matrix, and it follows 
that the lattice spanned by the <£(ci(Y, a)) is saturated, so it must equal the 
flow lattice of S. □ 

If S is a graph, then all its subcomplexes and relative complexes are 
torsion- free (equivalently, its incidence matrix is totally unimodular). There- 
fore, Theorems |6.1| and |6.2| give integral bases for the cut and flow lattices 
respectively. These are, up to sign, the integral bases constructed combina- 
torially in, e.g., [18j Chapter 14]. 

7. Groups and Lattices 

In this section, we define the critical, cocritical, and cutflow groups of a 
cell complex. We identify the relationships between these groups and to the 
discriminant groups of the cut and flow lattices. The case of a graph was 
studied in detail by Bacher, de la Harpe and Nagnibeda [2] and Biggs [I], 
and is presented concisely in [181 Chapter 14]. 

Throughout this section, let S be a cell complex of dimension d with n 
facets, and identify both C d (T,;Z) and C d (£;Z) with Z n . 

Definition 7.1. The critical group of S is 

K(Ti) := T(keTd d ^/imd d d* d ) = T(coker(im3 rf ^)). 

Here and henceforth, all kernels and images are taken over Z. 

Note that the second and third terms in the definition are equivalent 
because kerd d ^i is a summand of Cd_i(X;Z) as a free Z-module. This 
definition coincides with the usual definition of the critical group of a graph 
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in the case d = 1, and with the authors' previous definition in [15] in the 
case that £ is Q-acyclic in codimension one (when ker dd-i/ im d d d d is its 
own torsion summand). 

Definition 7.2. The outflow group of £ is Z n /(C(S)© 7(E)). 

Note that the cutflow group is finite because the cut and flow spaces are 



orthogonal complements in W 1 (Proposition 5.1), so in particular C T 
spans W 1 as a vector space. Observe also that the cutflow group does not 
decompose into separate cut and flow pieces; that is, it is not isomorphic to 
the group G = ((Flow d nZ n )/ker<9 d ) © ((Cut d nZ n )/im^), even when S is 
a graph. For example, if T, is the complete graph on three vertices, whose 
boundary map can be written as 







-1 




H 











^ o 1 


1 





then ker d = span{(l, 1, 1) T } and im<9 = span{(l, —1, 0) T , (1, 0, — 1) T }. So 
G is the trivial group, while Z n /(ker d d © im<^) = Ki(E) = Z 3 . 

In order to define the cocritical group of a cell complex, we first need to 
introduce the notion of acyclization. 

Definition 7.3. An acyclization of £ is a (d + l)-dimensional complex 
such that = £ and H d+ i(£l; Z) = H d (£l] Z) = 0. 

Algebraically, this construction corresponds to finding an integral basis 
for ker<9d(£) and declaring its elements to be the columns of d d +\ (^2) (so in 
particular |f2(d+i)| = /3d(E)). The definition of acyclization and equation ([!]) 
together imply that H d+1 (Q;7j) = 0; that is, d d+1 (£l) is surjective. 

Definition 7.4. The cocritical group K*(Y^) is 

#*(£) := C d+1 (Q;Z)/imd* d+1 d d+1 = cokerL^. 

It is not immediate that the group K*(J1) is independent of the choice 



of VL\ we will prove this independence as part of Theorem 7.7. For the mo- 
ment, it is at least clear that K*(Ti) is finite, since rank<9^ +1 = rankL^K = 
rankCd+i(f2; Z). In the special case of a graph, the cocritical group coin- 
cides with the discriminant group of the lattice generated by the columns of 
the "intersection matrix" defined by Kotani and Sunada [22J. (See also [5, 
Sections 2, 3].) 

Remark 7.5. As in |15j . one can define critical and cocritical groups in every 
dimension by ^(S) = T(C;(E;Z)/ imd i+1 d* +1 ) and^*(S) = T(C;(£; Z)/ imd*di 
If the cellular chain complexes of S and ^ are algebraically dual (for exam- 
ple, if £ and ^ are Poincare dual cell structures on a compact orientable 
(i-manifold), then K^) = K^iE) for all i. 

We now come to the main results of the second half of the paper: the 
critical and cocritical groups are isomorphic to the discriminant groups of the 



CUTS AND FLOWS OF CELL COMPLEXES 



23 



cut and flow lattices respectively, and the cutflow group mediates between 
the critical and cocritical groups, with an "error term" given by homology. 

Theorem 7.6. Let S be a cell complex of dimension d with n facets. Then 
there is a commutative diagram 



5- Z n /(C © J 7 ) 



■4> 







(14) 



K(T.) ^T(# d _i(E;Z)) 



> im d d j im d d d* d - 

in which all vertical maps are isomorphisms. In particular, K(Y,) = C^/C. 



Proof. Step 1: Construct the bottom row of (14). The inclusions \mddd*, C 



0. 



\m.dd C ker^-i give rise to the short exact sequence 

— > ixndd/ imdddj — > kerdd-i/ imddd d — > kerc^-i/ im9^ 

The first term is finite (because ranked = rank ddd d ) , so taking torsion 
summands yields the desired short exact sequence. 

Step 2: Construct the top row of ( fl4| ). Let r = ranked, let {v\, . . . ,v r } 
be an integral basis of C, and let V be the matrix with columns v\, . . . ,v r . 

the dual basis {i>*, ...,!)*} for C' consists of the columns 
V{V T V)~ 1 . Let ip be the orthogonal projection W 1 — > 

WV T = V{V T V)- l V T (see 



By Proposition 



2.2 



of the matrix W 
Cut(S), which is given by the matrix P 



Proposition 2.2). Then 



im ifi = colspace [v* 



[vi 



The i th column of P equals W times the i th column of V T . If we identify 
C s with 17 via the basis {v^, . . . ,v*}, then im(^) is just the column space 
of V T . So C/im^ = Z7colspace(F T ), which is a finite group because 
ranky = r. Since the matrices V and V T have the same invariant factors, 
we have 

Z r / colspace(y T ) T(Z"/ colspace(y)) = T(C d (Z;Z)/ imd* d ) = T(H d {Z;Z)). 

Meanwhile, im tp D C because PV = V . Since ker -0 = J~, we have (im i/))/C = 
(Z n /T)/C = Z n /(C © J 7 ). Therefore, the inclusions CCim^CC 1 give rise 



to the short exact sequence in the top row of (14). 



Step 3: Describe the vertical maps in (14). The maps a and j3 are each 



induced by dd in the following ways. First, the image of the cutflow group 
under d d is 

d d (Z n /(T®C))=d d (Z n /(kerd d ®imd* d )) = imd d /imd d d* d . 

On the other hand, dd acts injectively on the cutflow group (since the latter 
is a subquotient of Z n /ker dd). So the map labeled a is an isomorphism. 
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The cellular boundary map dd also gives rise to the map j3 : C^/C — > K (£), 
as we now explain. First, note that daC^ C imjR dd Q ker^c^-i- Second, 
observe that for every w G C and p G Cd_i(£;Z), we have (dw,p) = 
(w, d*p) G Z, by equation Q and the definition of dual lattice. Therefore, 
<9 d C C Cd_i(E;Z). It follows that d d maps C to (ker R nC d _i(S; Z) = 
ker^ hence defines a map /3: C/C — >• ker^ d d -i/ im^ ddC^. Since C/C is 
finite, the image of /3 is purely torsion, hence contained in K(E) . M oreover, 
ft is injective because (ker d d ) D C = J" n C = by Proposition f 



5.1 



Every element of M n can be written uniquely as c + / with c G Cut(E) 
and / G Flow(E). The map tp is orthogonal projection onto Cut(£), so 
d d (c + f) = ddC = <9d(V>(c + /))• Hence the left-hand square commutes. The 
map 7 is then uniquely defined by diagram-chasing. 

The snake lemma now implies that ker 7 = 0. Since the groups T{Hd-\ (E; Z)) 
and T(i? rf (£;Z)) are abstractly isomorphic by equation 0, in fact 7 must 
be an isomorphism and coker7 = as well. Applying the snake lemma 



again, we see that all the vertical maps in (14) are isomorphisms. □ 



Theorem 7.7. Let E be a cell complex of dimension d with n facets. Then 
there is a short exact sequence 

-»■ T(iT d _i(E;Z)) -»■ Z n /(C ®T) ^PjT ^ 0. (15) 

Moreover, #*(£) = J^jT. 

Proof. Let be an acyclization of £. By construction, the columns of the 
matrix A representing dd+i(Q) form an integral basis for T = ker <9 d . Again, 
the matrix Q = A(A T A)~ l A T represents orthogonal projection M. n — > Flow(E). 
The maximal minors of A have gcd 1 (because T is a summand of Z n , so 



the columns of A are part of an integral basis), so by Proposition 2.2 the 
columns of Q generate the lattice J-K Therefore, if we regard Q as a map 
of Z- modules, it defines a surjective homomorphism Z n — > J 7 ". This map 
fixes J- pointwise and its kernel is the saturation C := (C ® M) (1 Z". So we 
have short exact sequences — > C — ^ Z n / T — ^ / T — ^ and 

-> C/C -> Z n /(C 0. 



Since C is a summand of Z n by Proposition 2.1, we can identify C/C with 
T(F d (E;Z)) ^ T(^_x(S;Z)), which gives the short exact sequence jig. 

We will now show that J 7 " /J 7 S iT* (E). To see this, observe that d* d+1 ( T^) -- 
<9*j +1 (colspace((5)) = colspace(^4 T Q) = colspace(^4 T ) = imi9^ +1 = Cd+i(Q) 
(by the construction of an acyclization). In addition, keri9^ +1 is orthogonal 
to J 7 ", hence their intersection is zero. Therefore, <9^+i defines an isomor- 
phism J 7 ^ — y Cd + i(n). Moreover, the same map d* d+1 maps T = kerdd = 
imdd+i surjectively onto imd d+l dd+i- □ 

Corollary 7.8. If ff^_i(E;Z) is torsion-free, then the groups K(E), K*(E), 
C/C, J^jT, and Z n / '(C © J 7 ) are a// isomorphic to each other. 
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Corollary |7,8| includes the case that E is a graph, as studied by Bacher, 
de la Harpe and Nagnibeda [2] and Biggs [4] . It also includes the combina- 
torially important family of Cohen-Macaulay simplicial complexes, as well 
as cellulations of compact orientable manifolds. 

Example 7.9. Suppose that fld(E;Z) = Z and that fl^— i(E;Z) is torsion- 
free. Then the flow lattice is generated by a single element, and it follows 



from Corollary [7jj that K(E) = if*(E) = F^/F is a cyclic group. For 
instance, if E is homeomorphic to a cellular sphere or torus, then the critical 
group is cyclic of order equal to the number of facets. (The authors had 
previously proved this fact for simplicial spheres [151 Theorem 3.7], but this 
approach using the cocritical group makes the statement more general and 
the proof transparent.) 

Example 7.10. Let £ be the standard cellulation e° U e 1 U e 2 of the real 
projective plane, whose cellular chain complex is 

z z ^> z. 

Then C = imd% = 2Z, C» = ±Z, and K(E) = C*/C = Z 4 . Meanwhile, 
F = JH = F$/F = 2T*(E) = 0. The outflow group is Z 2 . Note that the 



rows of Theorem 7.6 are not split in this case. 



Example 7.11. Let a, b £ Z\{0}. Let £ be the cell complex whose cellular 
chain complex is 

z d2=[a % z ^ z. 

Topologically, £ consists of a vertex e°, a loop e 1 , and two facets of dimen- 
sion 2 attached along e 1 by maps of degrees a and 6. Then 

C t /C = Z T , Z 2 /(C(BF) = Z r/g , F i /F = Z r/g2 , 

where t = a 2 + b 2 and g = gcd(a, b). Note that r = T2(S) is the complexity 
of S (see equation Q) and that g = |JETi(E;Z)|. The short exact sequence 



of Theorem 7.7 is in general not split (for example, if a = 6 and b 



8. Enumeration 

For a connected graph, the cardinality of the critical group equals the 
number of spanning trees. In this section, we calculate the cardinalities of 
the various g roup inva riants of E. 



Examples 7.10 and 7.11 both indicate that K(T,) = C^/C should have 
cardinality equal to the complexity t(£). Indeed, in Theorem 4.2 of [15] . 
the authors proved that |if(E)| = t(£) whenever E has a cellular spanning 
tree T such that Hd-i(T; Z) = i^-i(E;Z) = (in particular, E must be 
not merely Q-acyclic, but actually Z-acyclic, in codimension one). Here we 
prove that this condition is actually not necessary: for any cell complex, 
the order of the critical group K(T,) equals the torsion-weighted complexity 
r(E). Our approach is to determine the size of the discriminant group C"/C 
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directly, then use the short exact sequences of Theorems 7\6 and 7/7 to 
calculate the sizes of the other groups. 



Theorem 8.1. Let X be a d- dimensional cell complex and let t = 
IT^x^Z))). Then 

|cVc| = |tf(£)|=T d (E), 

|z"/(ce^)|=T d (s)/t, 

\P/F\ = \K*(Z)\=T d (Z)/t 2 . 



and 



Proof. By Theorems 7.6 and 7.7 it is enough to prove that \C$/C\ = t^(S). 

Let R be a set of (cZ — l)-cells corresponding to a row basis for d (hence a 
vector space basis for Cut(S)); let 1Z be the lattice spanned by those rows 
(which is a full-rank integral sublattice of C); and let Y = (Srf_i\i2)US(^_ 2 ). 
The inclusions K C C C C d (E) give rise to a short exact sequence — > 
C/K -»■ # d (S,r ; Z) -> H d (Y,;Z) -> 0. Since C/ft is finite, the torsion 
summands form a short exact sequence (see Section [3]) . Taking cardinalities 
and using equation Q, we get 



\C/K\ 



|T(^ d (s,r ; z))| t d _!(s,r) 



(16) 



|T(i^(£;Z))| td_i(E) 

The inclusions ft C C C C s C ft* give |ft tt /^l = |^ tt /C tt | • |C tt /C| • |C/ft|. 
Moreover, TZ^/C^ = C/ft. By equation (16) and Binet-Cauchy, we have 



|C»/C| 



\0/TZ\ 



|^/c s l • |c/ft| 

t rf -i(S) 2 

t rf -i(s,r) 2 



t rf -i(S) 2 

td-i(s,r) 2 

det(^) 



l^/^l 



TCS d : |T| = 



t^l(S) 2 

;s,r)< 



det(<9p T ) 



By Proposition 3.2, the summand is nonzero if and only if T is a cellular 
spanning forest. In that case, the matrix 8r^ is the cellular boundary 
matrix of the relative complex (T,T), and its determinant is (up to sign) 
t<2-i(T, r), so by Proposition 3.4 we have 



\C*/C\ 



l(T) 5 



with the sums over all cellular spanning forests ICS. 



□ 



Dually, we can interpret the cardinality of the cocritical group as enumer- 
ating cellular spanning forests by relative torsion (co)homology, as follows. 



Theorem 8.2. Let £1 be an acyclization o/E. Then 

= J] \H d+1 (n, T; Z)| 2 = y \H d (Q, T; Z)| 2 

T 



T 
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with the sums over all cellular spanning forests ICS. 

Note that the groups H d+1 (£l, T; Z) and H d (tt,T;Z) are all finite, by 
definition of acyclization. 

Proof. Let d d +i = d d+ i(£l). Note that rank 9^+1 = f3 d (T,); abbreviate this 
number as b. By Binet-Cauchy, we have 

\K*CZ)\ = \detd* d+l d d+1 \= Yl ( detd *B) 2 

BCQ d : \B\=b 

where 8b denotes the submatrix of d with rows B. Letting T = Ti\B, we can 
regard ds as the cellular boundary map of the relative complex (Q, T), which 
consists of b cells in each of dimensions d and d+1. By Proposition |3.2[ the 
summand is nonzero if and only if T is a cellular spanning tree of Q/ d \ = £. 
(Note that the d+1, B,T,Q^ in the present context correspond respectively 



to the d, R, r, T of Proposition 3.2 ) For these summands, H d+1 (Q, T; Z) 



H d (p,, T; Z) is a finite group of order | det 8b\- O 

Then 



Remark 8.3. Let r*(£) = T ; z )| 2 > as in Theorem 

combining Theorems |8.1|and Theorems |8,2| gives 



8.2 



|C»/C| = |tf(£)|=T(£) =T*(S)-t 2 , 

I^V-^I = l^*( s )l = = -^(s)/t 2 , 

|Z n /(C0^)| =r(S)/t = r*(S)-t, 
highlighting the duality between the cut and flow lattices. 

9. Bounds on combinatorial invariants from lattice geometry 

Let n > 1 be an integer. The Hermite constant j n is defined as the 
maximum value of 

( mm <:r,*>)(|£V£|r 1/n (17) 

over all lattices C C W 1 . The Hermite constant arises both in the study 
of quadratic forms and in sphere packing; see |23[ Section 4]. It is known 
that 7„ is finite for every n, although the precise values are known only for 
1 < n < 8 and n = 24 [SJ. 

As observed by Kotani and Sunada [22J, if C = T is the flow lattice of a 
connected graph, then the shortest vector in T is the characteristic vector 



of a cycle of minimum length; therefore, the numerator in equation (17) is 
the girth of G. Meanwhile, \ J-^/J- \ is the number of spanning trees. We now 
generalize this theorem to cell complexes. 

Definition 9.1. Let E be a cell complex. The girth and the connectivity 
are defined as the cardinalities of, respectively the smallest circuit and the 
smallest cocircuit of the cellular matroid of S. 
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Theorem 9.2. Let S be a cell complex of dimension d with girth g and 
connectivity k, and top boundary map of rank r. Let b = rankij^-i(S; Z). 
Then 

fcr(S)- 1 /'- < 7r an d gT *(^)-^ b < 76 . 

Proof. Every nonzero vector of the cut lattice C contains a cocircuit in its 
support, so mux,, e c\ {0} ( x i x ) — ^- Likewise, every nonzero vector of the flow 
lattice J- of £ contains a circuit in its support, so min x gjr\ |Q|(x • x) > g. 
Meanwhile, |C/C| = r and iJ-^/J 7 ! = r* by Theorem 8.1 The desired- 
inequalities now follow from applying the definition of Hermite's constant 
to the cut and flow lattices respectively. □ 
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